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Abstract 


A brief  study  was  made  of  three  fallout  predic- 
tion models:  the  ENW  model  presented  by  Samuel 
Glasstone,  the  Miller  model  by  C.  F.  Miller,  and  the 
WSEG-10  model  by  George  Pugh  and  Robert  Galiano. 

Each  of  these  models  used  an  effective  wind  that  had 
constant  direction  and  speed.  A FORTRAN  computer 
code  of  the  Miller  model  was  prepared  by  the  author 
and  is  available  in  the  report. 

To  ascertain  the  effects  of  more  realistic  winds 
that  varied  direction  and  speed  with  altitude,  the 
author  developed  a model  that  utilized  an  altitude 
^ dependent  wind  as  well  as  a thin  stabilized  cloud, 

a log-normal  particle  size-activity  distribution,  a 
gaussian  distribution  of  activity  within  the  cloud, 
and  fall  time  equations  based  on  the  equations  of 
C.  N.  Davies.  This  model  was  prepared  as  a FORTRAN 
computer  code  by  the  author,  and  the  code  is  included 
in  the  report. 

The  two  most  significant  results  of  the  variable 
wind  model  are  the  asymmetric  pattern  produced  on  the 
ground  and  the  non-linear  centerline  of  that  pattern. 
The  model  allows  the  user  to  introduce  his  own  dis- 
cription  of  the  physical  processes  of  fallout  deposi- 
tion and  is  therefore  not  constrained  as  are  the 
f stylized  models  of  Glasstone  and  Miller. 

Vi 


A BRIEF  SURVEY  OF  FALLOUT  PREDICTION 
MODELS  AND  INTRODUCTION  OF  A FALLOUT 
PREDICTION  MODEL  UTILIZED  ALTITUDE 
DEPENDENT  WINDS 

I . Introduction 

The  most  significant  residual  effect  of  a land 
surface  nuclear  detonation  is  the  biological  hazard 
from  the  radioactive  debris  or  fallout  of  the 
explosion.  When  a surface  nuclear  detonation  takes 
place,  some  of  the  soil  nearest  the  explosion  is 
vaporized  by  the  intense  heat.  Mixed  with  the 
vaporized  soil  is  the  radioactive  residue  or  debris 
of  the  weapon.  This  debris  consists  of  the  remaining 
fissile  material  of  the  weapon  and  the  fission  pro- 
ducts of  that  portion  of  the  fissile  material  that 
had  fissioned.  This  mixture  is  carried  upward  within 
the  cloud  formed  by  the  detonation  and  at  once  beings 
to  cool  and  condense  into  particles  that  may  be  as 
small  as  a few  microns  (10~^  m) . Each  of  these 
particles  will  carry  some  of  the  radioactive  debris 
distributed  within  the  mass  of  the  particle  or 
condensed  on  its  surface,  and  each  will  eventually  be 
deposited  on  the  ground.  The  time  required  for  this 
deposition  will  vary  from  a few  minutes  for  the 
larger  particles  to  several  years  for  the  smaller 
particles . 
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The  deposition  of  this  radioactive  debris,  or 
fallout,  from  nuclear  explosions  has  been  a concern 
of  the  public  and  of  governmental  agencies  for  many 
years  because  of  the  immediate  health  hazards  and 
possible  long  range  genetic  damage  due  to  the  radia- 
tion. To  help  forecast  the  extent  and  the  degree  of 
contamination  of  the  earth's  surface  by  such  fallout, 
several  models  of  the  deposition  patterns  that  would 
predict  the  location  of  the  fallout  and  the  radiation 
intensity  were  prodcued  in  the  late  1950' s and  early 
1960's.  These  models  were  idealized  approximations 
to  experimental  data  gathered  from  American  weapons 
tests  that  included  the  effects  of  the  weapon's  total 
and  fission  yield,  and  the  mean  wind  speed.  The 
Miller  model  introduced  the  effects  of  fractionation, 
whereas  the  other  models  postulated  that  the  activity 
of  a fallout  particle  was  proportional  to  its  volume. 
The  factors  for  which  these  models  could  not  account 
include  variations  in  weather  conditions  such  as 
humidity  and  precipitation,  variations  in  terrain  and 
in  soil  content,  and  variation  of  wind  speed  and  wind 
direction  with  altitude. 

The  purpose  of  this  thesis  was  to  evaluate  three 
of  the  early  models  and  relax  one  of  their  constraints 
by  including  a variable  wind.  Local  fallout  will  be 
considered  as  that  fallout  deposited  within  1500  kilo- 
meters of  the  ground-zero  of  the  burst.  Such  factors 
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as  fractionation,  entrained  debris,  and  neutron  acti- 
vation of  the  soil  and  debris  will  not  be  addressed. 

The  author  will  briefly  examine  the  fallout  model 
presented  by  Samuel  Glasstone  (1963)  in  "The  Effects 
of  Nuclear  Weapons"  (ENW),  the  Miller  model  by  Carl 
F.  Miller  (1963),  and  the  WSEG-10  model  by  George  Pugh 
and  Robert  Galiano  (1959)  and  then  discuss  the 
developement  of  a model  produced  here  that  allows  for 
a horizontal  wind  that  varies  in  direction  and  speed 
with  altitude. 
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II.  The  ENW  Model 


This  model  was  first  presented  in  1957  (Ref  1) 
and  again  in  1963  (Ref  2)  after  the  incorporation  of 
new  and  more  extensive  information.  It  was  intended 
for  use  by  the  widest  possible  range  of  readers  and 
thus  offers  the  advantage  that  a high  level  of  tech- 
nical expertise  is  not  needed  for  its  use.  This 
model  is  presumably  an  empirical  fit  to  experimental 
data . 

Figure  1 (Ref  2:449)  displays  a typical  fallout 
pattern  as  presented  in  Reference  2.  This  pattern 
represents  the  unit-time  reference  dose-rate 
(Roentgens /hour  at  H+1  hour).  Reference  2 explains 
how  to  scale  this  pattern  for  other  yields  and  wind 
speeds.  The  simplicity  of  the  scaling  operations  for 
windspeed  and  weapon  yield  make  this  model  ideal  for 
use  by  a field  commander  with  limited  technical  assis- 
tance who  must  have  some  estimate  of  the  extent  and 
degree  of  fallout  contamination  so  that  he  can  limit 
the  exposure  of  his  men.  Several  charts  and  tables 
are  presented  in  Reference  2 that  further  enhance  the 
value  of  this  model.  These  include  protective  factors 
for  different  structures,  accumulated  absorbed  dose 
as  a function  of  time  exposed,  and  absorption  or 
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or  attenuation  coefficients  for  gamma  radiation  in 
various  materials. 

This  model  ignores  neutron  induced  activity,  stem 
fallout,  and  throwout  within  a blast-damage  circle 
about  ground  zero.  It  assumes  that  the  wind  will 
have  a constant  velocity  and  direction  and  will  re- 
main so  for  the  lifetime  of  the  deposition  process. 

As  previously  stated,  the  patterns  produced  by  this 
model  are  idealized,  therefore  the  reader  would  not 
expect  them  to  closely  approximate  the  experimental 
data  from  any  particular  burst.  This  author  has 
found  no  experimentally  obtained  fallout  pattern 
which  is  accurately  predicted  by  this  model.  Prepar- 
ing a computer  code  based  on  this  model  would  not  be 
difficult  if  some  accurate  figures  of  the  examples 
in  Reference  2 could  be  obtained.  If  this  data 
could  be  found  or  extracted,  then  the  resultant  com- 
puter code  would  be  quite  fast.  If  not,  then  this 
model  would  still  be  very  useful  in  a "handbook" 
status . 
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The  Miller  Model 


Carl  F.  Miller  first  published  the  results  of  his 
fallout  modeling  efforts  in  1963  (Ref  3).  The  Miller 
model  is  an  empirical  fit  to  experimental  data  that, 
like  Glasstone's  model,  utilizes  a constant  speed 
and  constant  direction  wind.  In  his  model  Miller 
included  the  results  of  investigations  into  the 
thermodynamics  of  fallout  particle  formation,  fire- 
ball behavior,  fractionation  of  fission  products, 
the  effect  of  wind  shear,  and  the  biological  effect 
of  ionizing  radiation. 

The  fallout  pattern  or  footprint  resembles  the 
shadow  of  a mushroom  shaped  cloud  characteristic  of 
.luclear  explosions.  This  cloud  is  described  by  the 
Miller  model  as  a truncated,  inverted  exponential 
cone  (stem)  topped  by  an  oblate  spheroid  (cloud). 

Of  the  models  surveyed,  the  Miller  model  was  unique 
in  two  respects.  First,  it  attempted  to  model  stem 
fallout,  and  secondly,  it  constructed  the  predicted 
fallout  pattern  around  several  characteristic  points 
of  location  and  dose  rate.  The  radiation  intensity 
or  unit-time  reference  dose-rate  varies  as  an  inverse 
exponential  with  distance  away  from  and  along  the 
pattern  centerline.  The  scaling  procedures  for 
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wlndspeed  and  weapon  yield  are  not  as  simple  here  as 
in  the  ENW  model,  but  then  the  scaling  procedures 
for  Miller's  model  are  of  a different  type  than  are 
those  of  the  ENW  model.  The  complexity  of  the 
Miller  model  does  not  permit  the  linear  scaling  used 
by  ENW,  however  the  Miller  model  does  more  adequately 
describe  the  physical  process  of  fallout  particle 
formation  and  deposition. 

Because  this  model  presents  an  idealized  pattern 
and  isointensity  contours  one  should  not  expect  it 
to  accurately  predict  the  fallout  deposition  from  a 
particular  burst,  but  to  present  a generalized 
approximation  to  many  actual  patterns.  The  pattern 
predicted  by  this  model  has  a shape  similar  to  that 
predicted  by  ENW  but  is  considerably  shorter  and 
thus  encloses  less  land  area  within  any  isointensity 
contour . 

The  thoroughness  and  completeness  of  this  model 
is  exemplified  by  the  variety  of  ways  available  to 
extract  data  from  it.  One  can  compute  the  unit-time 
reference  dose-rate  at  a point,  describe  an  isointen- 
sity contour  as  a function  of  two-directional 
displacement,  and  determine  the  area  within  a given 
contour.  This  author  chose  to  prepare  a FORTRAN  IV 
computer  code  of  this  model  based  on  the  adaptation 
of  it  presented  in  Reference  7.  Figure  2 presents 
the  output  from  this  code  for  a 1 MT , 100  percent 
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oi 

** 

in 

44 

•t 

44 

44 

44 

44 

ro 

44 

44 

44 

44 

44 

o 

(M 

44 

44 

44 

44 

t4 

• 

44 

44 

44 

44 

44 

44 

(VJ 

«• 

44 

44 

44 

44 

44 

44 

4^ 

to 

(T 

44 

44 

44 

44 

44 

(VJ 

OD 

44 

44 

44 

44 

44 

44 

44 

44 

44 

N. 

44 

44 

44 

44 

44 

44 

44 

44 

44 

«x> 

▼4 

44 

44 

44 

44 

44 

44 

t4 

in 

4^ 

44 

44 

44 

44 

44 

44 

44 

44 

44 

44 

44 

44 

44 

44 

44 

44 

44 

K). 

44 

44 

44 

44 

44 

44 

44 

44 

4-4 

o 

(Vi 

44 

44 

44 

44 

4H 

44 

44 

44 

44 

44 

44 

• 

▼4 

44 

44 

44 

44 

44 

44 

44 

44 

44 

CVJ 

44 

44 

44 

44 

44 

44 

44 

44 

44 

44 

44 

cr 

(T 

4^ 

44 

44 

44 

44 

4^ 

44 

44 

44 

44 

44 

«o 

44 

44 

44 

44 

44 

44 

44 

44 

44 

y* 

4< 

44 

44 

44 

44 

44 

44 

»4 

t4 

44 

44 

«0 

t4 

44 

44 

44 

44 

44 

44 

44 

44 

44 

44 

44 

44 

in 

44 

44 

44 

44 

44 

44 

44 

44 

44 

44 

44 

44 

j- 

t4 

4-i 

44 

4< 

44 

41 

44 

4H 

44 

44 

44 

44 

44 

fO 

44 

44 

44 

44 

44 

44 

44 

44 

44 

44 

44 

o 

(VJ 

44 

44 

44 

44 

44 

44 

44 

44 

44 

■r* 

44 

• 

▼4 

44 

44 

44 

44 

(VI 

44 

44 

44 

44 

t4 

CVi 

▼4 

44 

44 

4^ 

44 

(VJ 

(VI 

(VJ 

44 

44 

44 

44 

44 

uv 

O' 

44 

44 

44 

(VI 

(VJ 

(VJ 

44 

44 

44 

44 

▼4 

44 

44 

44 

44 

(VI 

tJ 

(VI 

44 

44 

44 

44 

44 

44 

44 

44 

(V.f 

(VJ 

(VI 

(VJ 

(VI 

44 

44 

44 

44 

vO 

44 

44 

44 

tJ 

(VJ 

to 

(VJ 

(v: 

44 

44 

44 

44 

ir»- 

• ^ 

44 

44 

44 

(VJ 

CVI 

CsJ 

(^ 

CV» 

44 

44 

44 

44 

▼4 

44 

44 

44 

(Vi 

(VI 

(VI 

(V, 

tJ 

44 

44 

44 

44 

44 

fr> 

44 

▼4 

44 

(VJ 

(VJ 

(VI 

CVI 

(VI 

44 

44 

44 

44 

4# 

o 

(VJ 

▼4 

44 

to 

(VJ 

(Vt 

ro 

twi 

(VJ 

(VJ 

44 

44 

44 

44 

• 

44 

44 

to 

(VI 

r^j 

ro 

ro 

n- 

(Vi 

44 

44 

44 

44 

CVi 

♦ 

44 

44 

(VJ 

tJ 

f^} 

ro 

ro 

M 

(VI 

44 

44 

44 

44 

IT 

»4 

44 

4< 

tJ 

M 

K' 

ro 

(VI 

nj 

44 

44 

44 

*30 

44 

r-4 

to 

(VJ 

Kj 

TO 

tJ 

(VJ 

44 

44 

44 

w4 

▼4 

44 

44 

(VI 

ro 

ro 

r^ 

fO 

ro 

(VJ 

44 

44 

44 

«£ 

44 

44 

44 

(VI 

ro 

ro 

ro 

ro 

ro 

(VJ 

44 

44 

44 

r4 

in 

44 

44 

44 

(VI 

ro 

ro 

ro 

ro 

ro 

(VJ 

44 

44 

44 

44 

44 

CVJ 

r»/ 

j- 

!•> 

ro 

(VJ 

44 

44 

t4 

ro 

▼4 

44 

(VJ 

(VI 

ro 

ro 

r^) 

ro 

(VJ 

(Vi 

44 

44 

44 

UJ 

(V. 

▼4 

44 

0'^ 

(VJ 

t. 

>;i 

•i 

h • 

i\j 

(VJ 

44 

44 

4^ 

• 

▼4 

44 

(VJ 

(VJ 

ro 

a 

ro 

(Vi 

(VJ 

44 

44 

44 

(VI 

* 

r4 

▼4 

44 

(Vi 

(VJ 

r'l 

■t 

j- 

ro 

(Vi 

(VI 

t4 

44 

44 

O' 

t4 

44 

(VJ 

(VI 

r»j 

■t 

-tf- 

r> 

(VJ 

(VI 

44 

44 

44 

cO 

44 

44 

ev. 

t'i 

ro 

-T 

^ 

ro 

(U 

4< 

44 

t4 

N. 

44 

4» 

'VJ 

ro 

r^j 

j- 

ro 

ro 

(VJ 

▼4 

44 

44 

vC 

44 

44 

(VJ 

tr> 

r^ 

t 

■t 

ro 

r*} 

t* 

44 

44 

44 

in 

«>4 

44 

44 

to 

ro 

ro 

in 

ro 

to 

44 

44 

44 

4H 

44 

(VJ 

fO 

-t 

uv 

j- 

.4 

t 

1 1 

44 

44 

44 

K 

44 

O' 

ro 

j- 

L"' 

^4* 

ro 

(V! 

44 

44 

44 

o 

(Vi 

44 

44 

(VJ 

ro 

.4 

tr 

if 

t; 

n 

44 

44 

44 

• 

4< 

44 

OJ 

ro 

j- 

b 

ro 

<v- 

44 

44 

44 

(VJ 

» 

44 

44 

tJ 

ro 

j- 

ir 

r^- 

(VI 

44 

44 

K' 

O'- 

• ^4 

44 

44 

tJ 

ro 

j- 

-4 

IT 

fo 

tg 

44 

44 

4t 

44 

44 

44 

(VI 

ro 

-J- 

ir. 

ro 

(Vi 

44 

44 

44 

44 

44 

44 

(VJ 

ro 

J- 

in 

(VI 

44 

44 

44 

«o 

a 

lA 

in 

ro 

to 

ro 

44 

to 

44 

in 

44 

o 

(VI 

44 

in 

• 

rt 

«r 

• 

» 

O O 

o 

o 

o 

o 

o 

O 

O 

O 

o 

C3 

O 

a 

(T 

cy 

c» 

(vi  ce 

O 

ir 

(VJ 

AT? 

rr 

cv» 

o o 

«r*  (Vj| 

Kj  (VI 

(VI  (VI 

44 

44 

1 

1 

44 

44  (Vj 

fVJ 

Cl  ro 

• 

1 

1 

1 

1 

u 

•o  o 

I—*  o 

0)  o 

•H 

>><U 

e " o 

C BS  -U 

o w X 

•H  I 

CO  05  O 

CO  1-t  on 

•w  C O 
CM  3 CO  I 

4J  <u  o -vt 
CXi  u 
4)  -U  •> 

O tH  oc 
Ui  «)  X 
<U  C I ^ — 
Q.-H  05 

3 tH 

O O O 

O X ••'O 

tH  C/i  TD  CO 

C 

•>  (U  o 

H • W)4J 

S'-v  4) 

i-l  i-l  O Dtf 

-H  4)  OX 

•a  D.^-c^ 


cd  o 

fl3 

Pi 

X 2 

1 

Vi 

o 

O M 

CO  o 

M-l  (U 

• 

o 

0> 

•» 

CO 

c ^ 

u 

b; 

CO 

X 

o 

0)  S 

u 

u 

U — ' 

u 

0) 

Pi 

o 

CO  4-> 

CO 

o o 

Q-  ui 

o 

o o 

U "O 

44 

44 

C 3 

O X 

01 

o 

1 

1-1 

o 

l-l 

u (U 

c 

m 

tH  u 

0) 

o 

CO  CO 

b 

CO 

A 

O 'W 

01 

Pi 

a.  bi  u-i 

1 

X 

01  3 

0) 

Q CO 

Wi 

CM  Pi 

• 

CM 

01 

V4 

3 

00 

•<-* 

b* 

(saiTW)  3uauiaDBX<^s -pQ  pu-pwssoaQ 


9 


fission  yield  burst  with  a 15  mph  wind.  The  most 
striking  dissimilarity  between  Figure  2 and  Figure  1 
is  the  "bud"  on  the  upwind  end  of  the  footprint  in 
Figure  2.  This  "bud"  represents  the  contribution 
of  stem  fallout.  One  can  readily  see  that  this 
contribution,  though  a high  dose  rate,  may  be  of 
slight  concern  to  those  persons  within  its  perimeter 
because  they  would  probably  be  dead  from  other 
effects.  A complete  listing,  and  glossary  of  terms 
of  this  computer  code  is  available  to  the  reader  in 
Appendix  A. 

The  Miller  model  as  presented  in  Reference  7 has 
several  singularities  at  which  the  model  fails  or 
gives  unreliable  results  (Ref  7:16-18).  These 
singularities  are  identified  by  certain  combinations 
of  weapon  yield  and  windspeed.  The  program  listing 
in  Appendix  A includes  a test  for  these  singularities 
in  the  subroutine  CONST.  The  user  must  observe  the 
limits  on  wind  speed  and  yield  of 

0 < wind  speed  < 75  mph  and  IKT  < yield  < 5000  KT . 
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IV.  WSEG-10 


Much  of  the  effort  put  into  the  preparation  of 
fallout  modles  was  motivated  by  the  operational  needs 
of  field  commanders.  WSEG-10,  first  published  in 
1959,  was  an  attempt  to  provide  such  commanders  with 
a tool  with  which  they  could  quickly  estimate  the 
location  and  radiation  intensity  of  radioactive 
fallout.  The  authors  of  WSEG-10,  George  Pugh  and 
Robert  Galiano,  had  the  following  to  say  about  the 
situation : 


"Fallout  estimates  for  use  in  operational 
planning  have  usually  been  obtained  either  by 
use  of  stylized  patterns  or  by  detailed  machine 
calculations.  Stylized  patterns  are  too 
inflexible  and  too  unrealistic  to  answer  many 
questions  encountered  in  operations  research. 
Detailed  calculations  which  have  been  used  pre- 
viously are  laborious  and  costly,  and  unless 
meteorological  conditions  are  known  in  extreme 
detail  they  do  not  produce  accuracy  of  results 
commensurate  with  the  effort.  The  purpose  of 
this  memorandum  is  to  introduce  a simplified 
computational  model  which  is  more  directly  tied 
to  the  physics  of  fallout  than  the  stylized 
patterns,  so  that  changes  in  physical  knowledge 
or  assumptions  can  be  more  readily  incorporated." 
(Ref  8:1) 


WSEG-10  was  an  attempt  to  describe  the  physical 
nature  of  fallout  by  modeling  the  spatial  distribution 
of  activity  within  the  radioactive  cloud,  assuming  a 
uniform  particle  activity  that  varied  only  with 
particle  volume,  and  describing  the  fall  of  each 
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particle  as  a sphere  in  a viscous  medium.  This  model 
describes  the  activity  distribution  within  the  cloud 
as  a normal  distribution,  and  the  particle  size- 
activity  distribution  as  a log-normal  distribution, 
as  does  ENW,  with  the  mode  radius  equal  to  28  microns, 
the  average  radius  equal  to  44  microns,  and  a 
standard  deviation  of  about  0.69.  The  use  of  a con- 
stant wind  direction  and  windspeed  was  probably  due 
to  its  intended  use  in  operational  situations  with 
little  time  and  limited  technical  assistance  availa- 
ble. 

The  results  of  this  model  were  not  presented  in 
"footprint"  form  as  were  the  results  of  Miller  and 
Glasstone,  but  were  given  in  tables  that  offered 
dimensions  and  locations  of  unit-time  reference 
dose-rate  isointensity  contours  as  functions  of 
particular  weapon  yields,  windspeeds,  and  wind  shear. 
Such  a presentation  could  be  very  valuable  to  a 
field  commander  in  an  operational  environment. 

Table  I is  an  example  of  the  tabular  presentation  of 
WSEG-10  results. 

The  patterns  produced  by  this  model  are  consid- 
erably longer  and  wider  than  those  of  Glasstone, 
indicate  peak  unit-time  reference  does-rates  well 
below  those  of  Glasstone  or  Miller,  and  encompass 
much  more  surface  area  than  either  of  the  two 
previous  models. 
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Table  I 

A Tabular  Presentation  of  WSEG-10  Results 
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The  advantage  of  WSEG-10  is  that  it  attempts  to 


predict  fallout  deposition  based  on  the  physical 
properties  of  the  atmosphere  and  the  fallout  parti- 
cles and  on  a statistical  analysis  of  the  particle 
distribution.  It  is  not  an  empirical  fit  to  experi- 
mental data.  This  attempt  to  predict  or  model  the 
physical  process  of  fallout  deposition  led  this 
author  to  develop  a fallout  model  that  would 
describe  the  physical  processes  of  falling  and 
dislocation,  and  that  would  allow  the  use  of  a wind 
that  varies  direction  and  speed  with  altitude. 


( 
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V.  A Variable  Wind  Model 

The  previous  models  each  used  a wind  that  was 
constant  in  speed  and  direction.  Only  Glasstone 
addressed,  but  did  not  incorporate,  the  fact  that 
atmospheric  winds  may  vary  in  both  direction  and 
speed  with  altitude.  This  author  developed  a fall- 
out deposition  model  that  would  incorporate  the 
effects  of  an  altitude  dependent  wind. 

There  were  several  simplifying  assumptions  made 
to  limit  the  scope  of  this  problem  and  form  the 
basis  for  its  solution.  These  are  listed  below: 

1.  the  source  of  all  fallout  particles  is  a 
thin  pancake  cloud, 

2.  all  particles  are  solid  spheres, 

3.  the  activity-particle  size  distribution  is 
a log-normal  distribution, 

4.  the  distribution  of  activity  horizontally 
across  the  pancake  cloud  is  a gaussian  dis- 
tribution , 

5.  there  is  no  fractionation  of  fission  products ; 
i.e.  the  activity  is  volume  or  mass  distributed 
within  each  particle. 

The  activity-particle  size  distribution  presented 
by  Miller  and  Sartor  (Ref  10:69)  was  chosen  to 
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represent  the  actual  distribution.  It  was  a log- 
normal distribution  with  a mean  radius  of  105 
micronSj  a mode  radius  of  30  microns,  and  a standard 
deviation  of  about  1.1.  This  distribution  was  ' 

separated  into  97  distinct  groups.  Each  group  was 
selected  so  that  the  largest  member  of  each  group 
would  fall  10  percent  faster  than  the  largest  member 
of  the  group  just  below  it  in  size.  The  smallest  I 

average  particle  radius  used  was  13.78  microns  for 
group  number  97.  The  group  containing  the  largest 
particles  was  group  number  1.  The  mean  particle 
radius  for  each  group  and  the  activity  fraction 
contributed  by  each  group  is  given  in  Appendix  B. 

Each  particle  group  was  treated  as  a separate 
cloud;  the  fall  time  for  which  was  computed  using  a 
drag-coefficient  Reynolds-number  method  (Ref  9:4) 
that  employed  the  density  and  viscosity  of  air  as 
given  in  Reference  11,  a particle  of  specific  gravity 
2.6,  and  Davies  .polynomials.  This  method  and  the 
Davies  polynomials  are  given  below  where  p is  the 
particle  density  (g/cm’),  is  the  air  density  (g/cm’), 
d is  the  particle  diameter  (cm),  n is  the  dynamic 
viscosity  of  air  (g/cm-sec),  g is  the  acceleration 
due  to  gravity  (constant)  (cm/sec*),  is  the  drag 
coefficient,  R is  the  Reynolds  number,  is  the 
terminal  velocity  (cm/sec),  AH  is  the  distance  fallen 
(cm),  and  T is  the  time  elapsed  to  fall  AH(sec). 
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(1) 


R 


4gp^  dV3n* 

- 2 . 3363X10 ' X ( C^RO"  + 2.0154X10  x (C^R*)’ 
- e.WOSXKj’xCC^RO"  for  < 138 


(2) 


Ios^qR  = - 1.29536  + 0.986  logj^QC^R*  - 0. 046677 (logj^QC^R^)* 

+ 0.0011235(logj^QC^RO’  for  138  <C^R*  <4.7X10"  (3) 

= (Rn/P3d)(l  + 2*^^^^/dp^)  (4) 

T = AH/V^  (5) 


The  second  term  in  parentheses  in  the  equation  is 
a correction  factor  for  small  particles  at  high 
altitudes  where  d and  p are  in  microns  and  grams  per 
cubic  centimeter  respectively.  The  time-to-fall  from 
various  altitudes  for  several  particle  sizes,  as 
computed  by  this  method,  are  very  much  like  those 
given  by  Glasstone  (Ref  2:496),  except  for  the  smaller 
particles  where  the  times  given  by  Glass.tone  are  much 
longer.  Figures  3 and  4 allow  some  comparison  of 
computed  fall  times  by  the  reader. 

The  height  of  the  radioactive  center  of  the 
stabilized  cloud  is  a function  of  weapon  yield  and 
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Figure  3.  Times  of  fall  of  particles  of  different 
sizes  from  various  altitudes  and  percen- 
tages of  total  activity  carried 
(Glasstone ) . 


was  taken  from  WSEG-10  (Ref  8:24)  in  the  form  of  the 
following  equation: 


H = 13.411  + 1.859  InY  - 0.0625(lnY  + 2.42)  |lnY  + 2.42 1 (6) 


where  H is  in  kilometers  and  Y is  the  weapon  yield  in 
megatons.  This  equation  is  graphed  in  Figure  5 
(Ref  8:25)  with  H in  thousands  of  feet.  The  equation 
for  the  cloud  radius  was  extracted  from  the  Miller 
model  (Ref  3:14)  and  is  the  following: 


R 


14. 661 X' 


:y0 


.431 


(7) 


18 


Time  (Hours) 

LI  for  particles  of  different  rad 
Ltudes  (Variable  Wind  Model). 


20 


Figure  5.  Cloud  height  as  a function  of  weapon  yield  as 
per  Equation  (6). 


where  R is  in  kilometers  and  Y is  in  megatons.  The 
cloud  height  is  used  to  determine  the  altitude  from 
which  particles  begin  their  descent.  Multiples  or 
fractions  of  the  cloud  radius  as  given  by  equation 
(6)  are  used  as  the  standard  deviation  of  the  parti- 
cles' horizontal  distribution.  This  distribution 
is  assumed  to  be  a circular,  symmetrical  gaussian 
distribution,  and,  because  it  describes  the  distri- 
bution of  each  particle  size  group,  it  can  be  used 
to  describe  the  distribution  of  the  activity 
contributed  by  each  particle  size  group.  A compari- 
son of  this  gaussian  distribution  of  activity  with 
the  typical  distribution  is  given  in  WSEG-10 
(Ref  8:23)  and  shown  in  Figure  6.  To  partially 
account  for  this  typical  distribution,  the  standard 
deviation  of  each  particle  size  group  was  increased 
by  0.020  of  the  standard  deviation  of  group  1 over 
that  of  the  particle  group  just  larger.  This  had 
the  effect  of  describing  a cloud  where  the  larger 
particle  sizes  were  more  concentrated  near  the  center. 
This  activity  distribution  is  shown  in  Figure  7.  Note 
that  this  distribution  is  flatter  than  the  regular 
gaussian  distribution  and  more  nearly  approximates 
the  actual  distribution  shape  given  in  Figure  6.  The 
flattening  of  the  activity  distribution  partially 
compensates  for  the  effect  of  wind  shear  by  dispersing 
each  particle  group  by  an  arbitrary  amount. 
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The  residual  activity  available  for  use  in  fallout 
formation  was  estimated  by  ENW  (Ref  2:492)  to  be 
approximately  550  gamma-megacuries  per  kiloton  of 
yield  at  H+1  hour.  The  average  energy  of  these  gamma 
rays  was  computed  to  be  0.902  MeV  (Ref  12:630), 
whereas  ENW  used  an  average  energy  of  0.95  MeV.  The 
dose  rate  was  calculated  for  a point  in  the  air  three 
feet  above  a smooth  infinite  plane.  For  a particular 
point  on  the  ground,  the  unit-time  reference  dose- 
rate  was  determined  by  summing  the  contributions  of 
each  particle  group.  To  compute  these  contributions 
it  was  necessary  to  know  the  location  of  the  center 
of  each  particle  group  cloud  on  the  ground.  These 
locations  were  determined  by  calculating  the  effect 
of  wind  direction  and  speed  on  the  group.  Displace- 
ments in  two  dimensions  were  calculated  by  multiplying 
fall  time  through  a wind  layer  one  kilometer  thick  by 
the  average  wind  speed  and  direction  in  that  layer. 

For  convenience,  wind  data  was  divided  into  20  layers, 
each  one  kilometer  thick.  The  reader  will  notice 
that  according  to  equation  (6)  the  weapon  yield 
required  to  produce  a stabilized  cloud  at  an  altitude 
of  20  kilometers  is  well  in  excess  of  100  MT  and  is 
much  larger  than  that  of  any  weapon  available  today. 

The  attenuation  of  the  gamma  radiation  in  the  air 
was  not  an  unusual  problem  as  the  following  derivation 
will  show.  The  flux  of  gamma  rays  at  point  X in 
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Figure  8 can  be  found  by  integrating 


^ (FLUX) 


A 2irr  “^air® 
41r»+  h*)  ® 


(8) 


to  yield 


FLUX  = 2.18284  A ■ ^-rays 

where  A is  the  gamma  activity  ( y-rays / sec-m* ) at  area 
a,  yCm”*)  is  the  macroscopic  cross  section  of  0.9  MeV 
gamma  rays  in  air  at  STP  conditions,  r(m)  is  the  ground 
distance  from  the  area  a to  the  point  for  which  the 
dose  rate  is  calculated,  and  s(m)  is  the  slant  range 
from  that  same  point  to  area  a.  The  dose  rate  is 
represented  by 

D = FLUX  X E X a „ (10) 

sec  - g 

where  o(m*/g)  is  the  linear  absorption  coefficient  of 
air  for  gamma  rays  of  average  energy  E (MeV).  The 
final  form  of  the  equation  for  dose  rate  is 

D = 1. 147125 X1(T“  a R/Hr  (§)  1 Hr  (11) 
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VI.  Results 


The  footprint  or  deposition  pattern  of  the 
variable  wind  model  for  a 1 MT,  100  percent  fission 
yield  burst  with  a constant  15  mph  wind  is  shown  in 
Figure  9.  The  similarity  between  Figure  9 and 
Figures  1 and  2 is  readily  apparent.  The  only  signif- 
f cant  differences  between  these  three  patterns  is  the 
variation  in  maximum  dose  rates,  the  variation  in 
contaminated  land  areas,  and  the  node  in  the  extreme 
downwind  portion  of  the  pattern  in  Figure  9. 

The  variable  wind  model  does  not  calculate  stem 
fallout  and  thus  ignores  what  could  be  a very  high 
dose  rate  area  near  ground  zero.  The  variation  in 
contaminated  land  areas  may  be  attributed  to  differ- 
ences in  fall  time  calculations  and  activity  distribu- 
tion. The  node  mentioned  above  is  due  to  the  lateral 
separation  of  particle  groups  and  is  caused  by  the 
markedly  increased  fall  times  for  the  smaller 
particles.  The  separation  of  particle  groups  can  be 
more  easily  seen  in  Figure  10  which  is  the  footprint 
for  an  18  KT , 96  percent  fission  yield  burst  with 
altitude  dependent  winds  as  given  in  Appendix  B, 

Table  II.  The  detached  "hot  spots"  of  radioactivity 
are  due  to  the  lateral  separation  of  the  smaller 
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East-West  Displacement  (Kilometers) 
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particle  groups.  The  effect  is  enhanced  by  the  scale 
used  on  each  axis.  A finer  or  smaller  scale  would 
show  areas  or  "splotches"  of  activity  instead  of 
activity  at  only  a few  points. 

The  asymmetry  of  the  pattern  in  Figure  9 clearly 
demonstrates  an  effect  of  altitude  dependent  winds. 
Previous  models  would  have  predicted  a symmetrical 
pattern  in  the  same  general  direction.  An  example 
of  the  source  of  this  asymmetry  is  given  in  Figure  11. 
The  figure  represents  the  actual  pattern  centerline 
for  a 10.4  MT , 70.2  percent  fission  yield  burst  with 
altitude  dependent  winds  as  given  in  Appendix  B, 

Table  III.  The  solid  curved  line  represents  the 
locus  of  the  particle  size  groups  as  each  impacted 
the  ground.  The  dashed  line  is  a straight  line 
connecting  the  location  of  the  first  and  last  particle 
size  groups.  This  figure  indicates  that  the  effects 
of  varying  altitude  dependent  winds  are  more 
pronounced  for  the  larger  particles. 

The  width  of  these  footprints  can  be  manipulated 
until  it  approximates  that  of  experimentally  obtained 
patterns  or  that  of  other  models  by  varying  the 
standard  deviation  of  the  particle  size  groups  or 
adjusting  the  particle  size-activity  distribution. 

The  downwind  extent  of  these  contours  is  still  con- 
siderably greater  than  that  of  other  models.  A FORTRAN 
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East-West  Displacement  (Kilometers) 

Figure  11.  The  locus  of  particle  groups  on  the 
ground  for  a 10.4  MT , 70.2  percent 
fission  yield  burst  with  varying 
winds  (Variable  Wind  Model). 


computer  code  of  the  variable  wind  model  and  a 
glossary  of  terms  is  available  in  Appendix  B. 
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VII.  Conclusions  and 
Recommendations 

The  usefulness  of  the  variable  wind  model  is 
obvious.  But  its  complexity  may  make  it  too  awkward 
for  use  in  a handbook  or  "yardstick"  fashion  and 
thus  unsuitable  for  field  commanders  in  an  operational 
environment.  If  technical  assistance  and  computer 
facilities  are  available  to  field  commanders  and  time 
varying  winds  are  employed,  this  model  may  provide 
a much  more  accurate  prediction  of  fallout  deposition 
than  any  previous  model. 

The  value  of  this  model  lies  in  its  potential  for 
describing  the  fallout  deposition  process.  To  improve 
the  existing  model  more  effort  and  attention  should 
be  given  to  determining  accurate  fall  times  and 
terminal  velocities,  describing  a more  realistic 
activity-particle  size  distribution,  more  accurately 
describing  the  activity  distribution  within  the 
stabilized  cloud,  and  to  estimating  the  residual 
gamma  ray  activity  from  a nuclear  burst.  To  expand 
the  model,  research  into  the  use  of  a cloud  with  a 
finite  thickness,  winds  that  vary  with  time  as  well 
as  altitude,  and  fallout  particle  formation  below  the 
stabilized  cloud  should  be  conducted  and  incorporated. 
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The  ripples  and  splotches  observed  in  the  pattern 
of  Figures  9 and  10  could  be  eliminated  by  increasing 
the  number  of  particle  size  groups  or  by  distributing 
the  activity  of  two  adjacent  groups  on  the  ground  so 
that  the  separation  between  each  group  is  filled  with 
some  activity.  This  could  also  reduce  the  downwind 
extent  of  each  isointensity  contour. 

This  model  offers  a tool  for  the  study  of  fallout 
deposition  that  should  be  used  and  explored. 
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Appendix  A 


The  Miller  Model 


The  required  and  unformated  input  for  this  pro- 
gram is: 


1.  total  weapon  yield  (MT) 

2.  fractional  fission  yield 

3.  mean  wind  speed  (mph) 

4.  the  number  of  specific  points,  if  any, 
for  which  the  user  desires  to  know  the 
unit-time  reference  dose-rate 

5.  the  number  of  expansions  of  the  pattern 
desired 

6.  the  coordinates  of  those  points  for  which 
the  user  desires  to  know  the  dose  rate, 
if  any  (miles) 


The  expansion  mentioned  in  5.,  above,  is  simply  a 
reduction  in  scale  of  the  first  printout  by  a factor 
of  one-third. 
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Glossary  of  Terms 


AB, 

Al,  A2,  A3, 

A4, 

A5, 

A6,  A7  - 

dummy  variables  used 

to  identify 

spe< 

:if  ic 

dose  rat 

es  in  the  fallout 

pattern . 

AJ 

- line  value 

of 

the  V 

ertical  axis  on  the  deposition 

pattern  (mi 

le  s ) 

• 

DX 

- increment 

used 

to  d 

itermine 

the  downwind  distance 

at  which  a ' 

dose 

rate 

is  to  be 

calculated  (miles). 

DXX 

- line  values  o 

f the 

hor izont 

al  axis  on  the 

deposition 

patt( 

arn  ( i 

miles)  . 

DY 

- same  as  fo 

r DX 

but 

for  cross 

wind  distance  (miles). 

FF 

- fractional 

f is 

Sion 

yield  of 

the  weapon. 

I - the  total  dose  rate  at  a specific  point 
(R/Hr.  @1  Hr.  ) . 

IBAK  - the  number  of  DX  increments  upwind  to  enclose 
the  upwind  portion  of  the  pattern. 

IC  - the  dose  rate  due  to  cloud  fallout  only 
(R/Hr.  @ 1 Hr.  ) . 

lEXP  - the  number  of  pattern  enlargements  desired. 

IFLD  - the  fallout  pattern  as  a rectangular  grid. 

IS  - the  dose  rate  due  to  stem  fallout  only 
(R/Hr.  @ 1 Hr.). 

ITEST  - a dummy  variable  used  to  check  for  discon- 
tinuities within  the  model. 

IX,  lY  - initial  values  of  DX  and  DY  respectively 
(miles ) . 
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123,  16,  17,  19  - dose  rates  at  the  characteristic 

points  X2 , X6 , X7,  X9  respectively  (R/Hr.  @ 1 Hr.). 
N - the  number  of  specific  points  for  which  the  user 
desires  dose  rate  values. 

PX,  PY  - the  downwind  and  crosswind  points  at  which 
a dose  rate  is  calculated  (miles). 

V - the  average  wind  velocity  during  fallout  deposi- 

tion (miles  per  hour),  0<V  <75. 

W,  WY  - the  fission  yield  and  total  yield  of  the 
weapon  (megatons). 

X - the  upwind  or  downwind  distance  to  a specific 

point  for  which  the  user  desires  a dose  rate 
(miles) . 

XI  - the  maximum  upwind  extent  of  the  IR/Hr.  contour 

(miles)  on  the  pattern  centerline. 

X2  - the  center  of  the  stem  fallout  pattern  (miles) 
on  the  pattern  centerline. 

X6 , X7  - the  distance  to  the  most  upwind  and  most 
downwind  extent  of  the  high  radiation  intensity 
ridge  due  to  cloud  fallout  (miles). 

X9 , X9P  - the  maximum  downwind  extent  of  the  IR/Hr. 
contour  for  a 15  mph  wind  and  any  other  wind- 
speed  respectively  (miles). 

Y - the  crosswind  distance  to  the  specific  point  for 

which  the  user  desires  a dose  rate  (miles). 

YS  - the  stem  pattern  half-width  of  the  IR/Hr.  contour 
at  X2  (miles ) . 
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Y6  - the  cloud  pattern  half-width  of  the  IR/Hr. 
contour  at  X6  (miles). 

Y8,  Y8P  - the  cloud  pattern  maximum  half-width  of 
the  IR/Hr.  controu  for  a 15  mph  wind  and  any 
other  windspeed  respectively  (miles). 
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PROGFftH  MILLri^dNPUT, OUTPUT) 

OlMfUSION  IFLO(1?0,23) ,0XX(13) 

COMMCN/Bv,OKl/tf,W,IT£?T,WV,FF/PLO<2/»X,PY,IS,IC/aLOX3/X2,X<»,X5,X6, 
1 xr,X9P,VS,Y8P,Y6,I23,I6,I’,r9 

REAL  Ills, 10,123,16., 17, 19 
INTEGER  A9,41,a2,A3,A4,ft5,A6,«7 

DATA  A9,Al,42,A3,Ai.,A5,46,a7/lH  , IHl,  1M2,  H3,  IHL  , 1H5,  IHfi,  1H7/ 

1 REAO*,HY,FF,V,N,IEXO 
XrCNY.LT.0.0)  STOP 
PR1MT82,WY,FF, V,M, ISXP 
ir((V.LE.O.O).OR.(Y.GT.75.))  GO  TO  199 
N»MY*FF 
CALL  CONST 

lEllTEST.EO.O)  GO  TO  150 

IPdtXP.GT.O)  GO  TO  200 

00  100  J>1,N 

REAO*,X,Y 

PX«X 

PY»A9S(Y) 

CALL  FIELD 
I»IS»IC 

too  PRINTSO,X,Y,I,IS,IC 
60  TO  1 

ISO  PRINT86,ITEST 
60  TO  1 
199  PRINTS  1 

GO  TO  1 . 

ZOO  PRINTS'^ 

IY«Y«P/22.»1 
IX»(X9P-X2»YS) /1?0.*1 
IFIIY.GT.IY)  IY=IX 
XFdX.LT.IY)  IX  = IY 
OY»XY 
OX*IX 

I8AX=(X2-YS)/0X-1. 

205  00  340  K=l,23 
1IC»2«.-K 
PY»{F-1)»0Y 
00  340  L«l,12a 
PX»<l-l*I94K)*0X 
CALL  FIELD 
I>IS«IC 

IFd-1.)  309,310,301 

301  lFd-30.)  310,311,302 

302  IFd-lCfl.)  311,312,303 

303  IFd-300.)  312,313,304 

304  IFd'lOOO.)  313,314,305 

305  IFd>3C0C.)  314,315,306 

306  IFd-lCOOO.)  315,316,316 

309  IFLn(L,IK)=AB 
GO  TO  340 

310  lFLO(L,IK)sAl 
60  TO  340 

311  ZFL0(L,I<)sA2 
60  TO  340 

312  1FLO(L,IK)pA3 
GO  TO  340 
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313  irLO(LtIK)sA% 

GO  TO  3<i0 

314  IFL0(L>I<>»A5 
GO  TO  340 

315  IFLn(L,IK>sA6 
GO  TC  340 

316  IFV.n<L,IK)sA7 
340  CONTINUE 

00  3e0  J*l,13 

360  OXX(J>=OX*tI96K*(J-l)*10.) 

PRINT83»0X, OY, DXX 
00  375  J=1,Z3 
AJ»(23-J)*0Y 

375  P'«NT64,AJ,(IFL0<M,J)  ,M=1,1Z0) 

00  376  J=1«ZZ 

IM»Z3-J 

AJ«-J*OY 

376  P«lNTe4,AJ, (IFLD(M,IM),M=1,120) 

IEXP=IEXP-1 

IF(IEXP)  1,1,380 
380  OY*Z./3.*OY 
0X*2./3.«0X 
GO  TC  Z05 

00  FORMAT  (“  “/5X-AT  THE  POINT  X*  HILE5,  AnO  Ya  ~F5.1 

1 -MILES,  THE  RADIATION  INTENSITY  DUE  TO  FALLOUT’VSX 

2 "IS  "F9.1-  R/HR,  with  "rq.f*  R/HR  DUE  TO  STEM  FALLOUT  AND' 

3 F9.1-  R/HR  DUE  TO  CLOUD  'ALLOUT”/) 

•1  FORMAT  r “//lO  {”X")ZX-HIND  SPEED  O'JT  OF  LIMITS,  PROCEED  TO  - 
1 "NEXT  PR09LEM“/A///) 

02  FORMAT <~1“4X”INPIIT  DATAl  "3F10.3,»T4) 

03  FORMAT C“1“5X“H0RI2CNTAL  (DOWMWINII  SCALF  IS  “F6.Z“  TO  1.  MILFS" 

1 5X  "VERTICAL  (CROSSWIND)  SULE  13  "F6.Z"  TO  1,  MILES"/// 

2 5X,F6.1,1Z(4X,F6. l> /PX 1 Z( "I Z 3455789*") //) 

04  FORMATC*  “F6. 1 , IX , IZ  0 Al ) 

06  FORMATC"  ~5X"ITEFT=  "IZ,  5X~OI'5S3NTINUITY  IN  CONSTANTS  OUF  TO  “ 

1 “INPUT  DATA  OR  MODEL  LIMITS.  PROCEED  TO  NEXT  OR03L EM. "///) 

07  FORMAT {“  “4X"MAP  LESFNO  "//SX"!-  1 TO  30  R/H»"/5X~2-  30  TO  100" 

1 “ R/HR"/5X"3-  ICO  TO  300  R/MR"/5X"4-  300  TO  1000  R/HR"/5X 

2 “5-  1300  TO  3000  R/HP"/5X"5-  3000  TD  10,000  R/HR"/5X 

3 "7-  10,000  OR  MORE  R/HR"//) 

ew 
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S08RCUTINE  CONST 

COHMCN/Pl.O<l/Vf  H,  ITEST,MY,FF/3L3<3/<2,X<*,i(5,X6,X7,X9P,  YS,Y8P,y6t 
1 123,16,17,19 

REAL  123,16,17,19 
H05>W/9. 

V20SV/20. 

I23««.0  00.*ll./W0^)*‘0.‘"2»<l./V?0»**a.75 
16*372 0.*S9PT( NOS ) *(l .-( V-25.) ••2/250  0.) 
I7*5C00.*S9RT('<(O=)^(i./V20J**Q.U0 

19*15. /W 

X2sl.75+0.23^( V-20 .)  •W05**0.23 
Xl**35.*V20^HO5^^C.23 
YS*7.1*HO5**0.35* (l./»20)*^0.75 
X1*X2-YS 

X6*2f .•V20*WO5**0.20 

X5*XF-11.*W05^^C. 30*(1.*(V-20.*W3;**0.30) /SO.* (V-20.*W05**0.30) 

1  ••2/600.) 

X7*67.1*H05^^0.30^tf20 
Y6*2f .•«!./ V20)^*0  .90 ♦H05*^ 0.30 
V8aA5.*(l./V2C)^*0.56*MO5**0.32 
X9a5E2.*H05“C  .3Q^V2C 

Y6PaVe*ALOG10 ( 17) / (ALOGIO (I7)-flL0313 (19) ) 

XgP*(X9*ALOG10 (17) -X7^4LOG10 (I9M /(4LOG10 ( 17) -ALOGIO (19) ) 
PRlNT85,wy, V,H,FF,Xl,X2,X(»,X5,X6,X/,X9P,r5, Y6, Y6P,I23,I6,I7,r9 
1TEST«0 

IFC  (123.GT.1.) .4tlO.<I6.GT.l.) . 4S0. <17 .GT . i . ) .ANO. (I7.G7 .19) .4N0. 

1 <X<,.GT.X2) .4N0. (X6.GT.X5)  .4N0. ( X7 . GT . X6) .4N9. (X90.GT.X7) 

2 ,*NC.  (X2.GE.0.0) .AN0.( YS.LE.(X4-X2)  ) .AND. (Y8P.LE. (X9P-X7)) ) 

3 ITE'T=1 

89  FORMAT (“  "//5X"PP09LEM  0AT4“//5X-.(E4P0N  YIPLO  - “F7.3~  NT-/5X 
1 ~M1N0  SPEED  - “F(,.l“  MPH-/5X-FISSI0N  YIEUO  - ~F7.3“  NT" 

•2  /5X“FRACT1CNAL  fission  YIE.D  - "05 . 3 ///5X , 60 ("*  'MZ/SX 

3 "••  NOTE  •*  ALL  distances  (X*S  AND  Y'S)  ARE  IN  -tLES" 

6 ",  ANO  ALL  INTENSITIES  (I'S)  ARE  IN  R/HR  AT  1HR."//5X 

9 60<~*  ") //'■X"OPOSINENT  points  0^  HE  P4TTERN“//GX“X1=  “F6.1 

6 ,3X“X2=  "Ft .1,3X"X4=  "ce.l.IX-XS*  "-5.1,3y~X6*  “C6.1,3X 

7 "X7»  "F6.1,3X“X9P=  ~F6.1/13<“YS=  5. l,29X"Y6s  "FS.i.IX 

8 ~Y8P*  “F5.1 ,/18X“I23*  "F9. 1, 2SX"I6*  "F7.1,2X~I7a  ~F7.l,2X 

9 "19*  “F7.1///) 

RCTUPM 

END 


f 
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SUBROUTINE  FIELD 

C0MMCN/BL0<2/PX,PY,IS,IC/BL0K3/X2,XinX'i,F5,Xr,Xgp,YS,Y8P,Y6,I23, 
1 16,17,19 

REAL  IS,IC,I23,I6,I7,I9 
ZSaO.O 
IC»0.0 

IF(PX-X4)  502,502,355 

502  IF(PY-YS)  503,503,599 

503  1P(PX-X2*YS)  599,515,505 
505  IF(PX-X2)  515,510,510 

510  IS=I23** (l.-SQRT ( ( (PX-X2) / ( X4-X2) ) •♦2» (PY/YS) ••2) ) 

60  TO  555 

515  IS=I23»* (l.-SORT ( (ox-X2)**2+PY*»2) /YS) 

555  IF(PX-X5)  599,557,356 

556  lF<Py-X9?)  557,557,599 

557  IF(PY-Y3P)  560  ,5f 0,599 
560  IF(PX-X6)  570,590,550 

570  IC=I6**(1.-S0RT(  ( (X5-PX)/(X6-X5))»*2MPY/Y6)**2»  ) 

RETURN 

500  IF(Py-X7)  590,590,595 
590  A=CX7-PX)/(X7-X6» 

6>1.-A 

IC» (I6**A*r7**85 (l.-PY/( A*Y6»9»Y8P) ) 

RETURN 

595  1C=I7**(1.-SQRT(( 0PX-X7I/(X9P-XTJ)**2»(PY/Y8P)**2)) 

599  RETURN 
END 
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Appendix  B 


The  Variable  Wind  Model 
Table  II 

97  Group  Particle  Size-Activity  Distribu- 
tion from  a Log-Normal  Distribution  with 
a Mean  Radius  of  105  Microns  and  a 
Standard  Deviation  of  0.69. 


HE  A'! 
GROUP 

RAorus 

(HIC’ONS) 

ACTIVITY 

FRACTION 

CONTRIOIJ- 

TFO 

MEAN 

GROUP 

PAOII'S 

(MICRCNS) 

ACTIVITY 

FPACTICN 

CONTRI0IJ- 

TEO 

mfak 
GROUP 
RAO I US 
(MICRONS) 

activity 

FRACTION 

C0NTRI9J- 

TSO 

13  V.  22 

.001191950 

291.02 

.011251860 

63.33 

.015^48350 

12%.  99 

.001316200 

277.48 

.011701620 

60.39 

.015227430 

1213.66 

.001450900 

264.56 

.0121463ri 

57.58 

.014885160 

113O.09 

.001596170 

252.25 

.012534753 

54. 9C 

.014523320 

1135. 1*. 

.001752911 

240.51 

.013014313 

52.34 

.014143700 

1033.71 

.001921211 

229.32 

.013433303 

49.91 

.013748180 

1 jni*.68 

.OC’101S50 

218.65 

.313039790 

47.58 

.013339660 

957.92 

.302293161 

208.47 

.0142ilS‘=3 

45.37 

.012917090 

913.34 

.002501550 

198.77 

.014607570 

43.26 

.01245538  0 

870.94 

.002721390 

199.52 

.014965110 

41.25 

.012045480 

8’0.31 

.002955000 

160.70 

.015302550 

39.33 

.011599290 

791.57 

.00.320264" 

172.29 

• 01561 8460 

37.50 

.011145630 

754.83 

.303464530 

164.27 

.015910910 

35.75 

.010695490 

719.70 

.OO’TAOSIO 

156.63 

.0161794*^0 

34.09 

.013241490 

695.21 

. 0040  ■'1350 

149.34 

. 016**1 965  0 

32.50 

.009758370 

654.27 

.004336750 

142.39 

.fll66331°0 

30.99 

.009337760 

623.82 

.00'4656290 

135.76 

.016817930 

29.55 

.009891190 

594.79 

.004990010 

129.44 

.016972930 

28.17 

.008450110 

567.11 

.005-»  37630 

123.42 

.017097050 

26.86 

.008015350 

543.72 

.005599750 

117.68 

.017189693 

25.61 

.007533660 

515.56 

.0060’290Q 

112.20 

.017250933 

24.42 

.007172550 

491.55 

.006459480 

106.96 

.017279530 

23.28 

.006763340 

469.69 

.006957790 

102.00 

.017275930 

22.20 

.005370150 

446.88 

.007267300 

97.25 

.017239750 

21.17 

.005936350 

4’6.08 

.007686230 

92.73 

.017171490 

20  .18 

.003615120 

435.25 

.009114340 

88.41 

.017071430 

19.24 

.005257050 

397.35 

.009550270 

84.30 

.016940140 

18.35 

.004912580 

369.32 

.008992770 

80.37 

.016778340 

17.49 

.004532070 

332.13 

.309443353 

76.63 

.0 16586930 

16.68 

.004265790 

335.75 

.009891670 

73.07 

.016366960 

15.90 

.00395339  J 

323.12 

.010345140 

69.67 

.016119530 

15.16 

.003676460 

315.22 

.310799110 

66.43 

.015846280 

14.46 

.003403470 

13.78 

.003144850 

y 
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Table  III 


Wind  Data  For  Figure  9 


Altitude 

(km) 

Wind  Direction 
(degrees ) 

Wind  Speed 
(km/hr) 

0-1 

90 

8.0 

1-2 

90 

25.7 

2-3 

90 

25.7 

3-4 

95 

27.4 

4-5 

115 

27.4 

5-6 

152 

22.5 

6-7 

152 

22.5 

7-8 

170 

24.1 

8-9 

170 

24.1 

9-10 

220 

32.2 

10  - 11 

220 

32.2 

11  - 12 

230 

27.4 

12  - 13 

230 

27.4 

13  - 14 

230 

27.4 

15  - 16 

2?0 

22.5 

16  - 17 

220 

22.5 

17  - 18 

220 

22.5 

18  - 19 

220 

22.5 

19  - 20 

220 

22.5 
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Table  IV 


Wind  Data  for  Figure  10 


Altitude 

(km) 

Wind  Direction 
(degrees ) 

Wind  Speed 
(kph) 

0-1 

66 

38.6 

1-2 

60 

37.0 

2-3 

79 

11.3 

3-4 

149 

4.8 

4-5 

134 

24.1 

5-6 

100 

29.0 

6-7 

100 

29.0 

7-8 

62 

16.1 

8-9 

184 

22.5 

9-10 

270 

27.4 

10  - 11 

270 

27.4 

11  - 12 

220 

59.5 

12  - 13 

220 

59.5 

13  - 14 

290 

56.3 

14  - 15 

290 

56.3 

310 

62.8 

230 

11.3 

17  - 18 

230 

11.3 

18  - 19 

260 

27.4 

19  - 20 

260 

27.4 
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Glossary  of  Terms 


A,  A2  - the  change  of  altitude  and  the  average  alti- 
tude respectively  through  which  a particle  falls 
(kilometers ) 

ACT  - the  total  residual  gamma— ray  activity  from  the 
burst  ( gamma ' s / sec ) 

AC,  AE , AF , AM,  AN  - variables  used  to  reduce  computer 
time  requirements  by  postulating  a linear  pattern 
centerl ine 

BA,  B1 , B2 , B3,  B4 , B5  - dummy  variables  used  to 
identify  specific  dose  rates  in  the  fallout 
pattern 

CDR  - the  dimensionless  variable  C^R*  used  to  compute 
a Reynolds  number  and  hence  a terminal  velocity 
for  each  particle  size 

CR  - the  stabilized  cloud  radius  as  given  by  Miller 
(kilometers ) 

CZ  - the  stabilized  cloud  height  as  given  by  WSEG-10 
(kilometers ) 

DB  - a dummy  variable  used  in  fall  time  calculations 

DH  - the  change  of  altitude  that  a particle  experi- 
ences as  it  falls.  The  same  as  A.  (kilometers) 

DR  - the  dose  rate  at  a particle  point  PX,  PY 
(R/HR  @ 1 HR) 

DT  - the  time  required  for  a particle  to  fall  a 
^ distance  A or  DH  (hours) 
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DX,  DY  - the  horizontal  displacements  either  east- 
west  or  north-south  experienced  by  a particle 
group  as  it  falls  through  any  wind  layer,  also 
the  DX  and  DY  increments  used  in  the  fallout 
pattern  determination  (kilometers  or  miles) 

DXL,  DYL  - the  number  of  DX  and  DY  increments, 

respectively,  below  or  behind  ground  zero  that 
would  allow  the  minimum  upwind  and  crosswind 
extent  of  the  fallout  pattern  to  be  included 
in  the  fallout 

DXX  - line  values  of  the  horizontal  axis  on  the 

V. 

deposition  pattern  (kilometers  or  miles) 

EX  - the  exponent  used  in  determining  each  particle 
group's  contribution  to  the  total  dose  rate 
FF  - the  fractional  fission  yield  of  the  weapon 
GDL  - the  ground  level  altitude  of  the  burst 
(kilometers ) 

H - the  same  as  AZ  (kilometers) 

I - a dummy  variable  used  to  determine  AZ 
lEXP  - the  number  of  pattern  enlargements  desired  by 
the  user 

IXY  - a single  line  of  the  pattern  for  a specific 
value  of  PY 

NP  - the  number  of  points  for  which  the  user  desires 
a dose  rate 

NN  - a dummy  variable  used  in  the  fall  time  determin- 
ations 
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PIR  - a conversion  factor  used  to  convert  degrees 
to  radians  ( radians /degree ) 

PR  - an  array  used  to  store  all  pertinent  information 
about  a particle  group:  average  particle  size 
(microns)  of  the  group,  fraction  of  total  activity 
contributed  by  this  group,  the  standard  deviation 
of  the  spatial  distribution  of  a group  (kilometers), 
the  time  to  fall  from  the  stabilized  cloud  to  the 
ground  (hours),  the  east-west  displacement  of 
the  group  cloud  center  (kilometers  or  miles),  and 
the  north-south  displacement  of  the  group  cloud 
center  (kilometers  or  miles) 

PX,  PY  - the  coordinate  of  a point  within  the  pattern 
for  which  a dose  rate  is  to  be  calculated 
(kilometers  or  miles) 

R - the  particle  radius  for  which  a fall-time  is 
computed  (microns) 

RE  - the  Reynolds  number  cal^'ulated  from  CDR  and  used 
to  compute  the  terminal  velocity 
SIGMA  - the  value  of  the  standard  deviation  of  the 
spatial  distribution  for  the  largest  group 
(kilometers ) 

SUM  - a dummy  variable  used  to  sum  the  dose  rate 
contributions  by  each  particle  group 
u - the  dynamic  viscosity  of  air  at  any  altitude 
(kg/m-sec) 

I 
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VT  - the  average  terminal  velocity  of  a particle  as 
it  falls  through  DH 

WIND  - an  array  that  stores  all  the  wind  and  altitude 
data:  the  average  true  direction  of  the  wind 
(360°  is  true  north)  that  is  later  converted  to 
machine  direction  (degrees)  for  each  kilometer 
of  altitude,  the  average  windspeed  (kph  or  mph) 
for  each  kilometer  of  altitude,  and  the  mid- 
point density  of  the  air  for  each  kilometer  of 
altitude  (kg/m*) 

YH  - a dummy  variable  used  to  limit  the  computer 
time  required  by  the  program 
YLD  - the  total  weapon  yield  (megatons) 

YLDO  - a dummy  variable  used  to  limit  the  computer 
time  required  by  the  program 
YND  - the  same  as  for  AC,  AE , AF,  AM,  and  AN 
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A FORTRAN  Computer  Code  for  the  Variable  Wind  Model 


The  required  and  unformated  Input  for  this  pro- 
gram Is: 


1.  the  mid-range  particle  radius  of  each 
particle  group  In  microns  and  the  fraction 
of  the  total  residual  activity  carried  by 
each  group 

2.  the  true  wind  direction  and  wind  speed  of 
each  1 km  thick  layer  of  atmosphere  up  to 
an  altitude  of  20  km 

3.  total  weapon  yield  (MT) 

4.  fractional  fission  yield 

5.  ground  level  of  the  area  contaminated  by  the 
fallout  (KM),  burst  Is  at  sea  level 

6.  the  number  of  specific  points  at  which  the 
user  desires  to  know  the  unit-time  reference 
dose  rate  (If  any) 

7.  the  number  of  pattern  expansions  desired 

8.  the  coordinates  of  those  points  from  4. 


The  expansion  mentioned  In  5.  above  uses  a scale  down 
factor  of  one-half. 
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WOG'»*M  NORHSL  (INPUT,  OUTPUT) 

COMMCN/9LOK1/WINO<2C , 3) , OT , OX, OY/3LOK 2/PR ( 97,6) , OR, ACT 
THE  PATTERN  IS  A “AXXHUH  OP  100  SPACES  WIOE 
OIHEHSION  IXY(IOC) ,0XX<11) 

INTEGER  PA, 81, 92, 93, 84, 95 

DATA  9A, 31, 32, 33, 84, 35  / IH  , IHl, 1H2, 1H3, 1H4, 1H5  / 

DATA  STATEHENT  LOAOS  HIO-POINT  AIR  DENSITY  -OR  EACH 
OF  ALTITUDE 

DATA  HIND<1,3) ,WTH0(2,3) ,HIND(3,3) ,NINO(4,3) ,WIN0(5,3) ,HINn(5,3) 

1 HIND(7,3> ,UIMO(8, 3) ,HIM0 (9,3)  ,HINO(ia, 3) ,HIN0(11,3) , 

2 HIN0(12,3) ,HINO(13,3) ,WIM9(14, 3) ,HIND(15,3) ,WIN0(15,3) , 

3 )(IN0(17,3)  ,«INO(  18,3)  , wind  (13,3)  , WIND  (20, 3) 

4 /I  .167  30,1  .05810, 0.95695,0.8  5343, 3 .7  7 7 04, 0 .6  97>,7, 0 . 62431 , 

5 C.  55719,0.  4957  6,  0 .4395^  , 0 . 383  57 , 0 . 33743 , 0. 28  53  3, 0 . 2'*646, 

6 0.21066,0.18036,0.15391,0.13157,0 .11248,0.09526/ 

READ*,  ((PR(N,H),M  = 1,2) ,N  = 1,97) 

101  READ*,  ((MINO(N,>(),H=l,2)  ,N  = 1,20) 

PRINTS 0, (((PR( N, H) ,4  = 1,3)  ,(HINO(N,L)  ,L=1,3))  ,N  = 1,20), 

1 ((PR(N,M),M=1,3) ,N=21,97) 

CONVERT  TRUE  HIND  DIRECTION  TO  HACRINE  DIRECTION 
DO  102  J=l,20 
NINO ( J,  1) =450. -KINO (J,l) 

IF(UINO(J, 1) .GT.360.)  WIND (J,l) ^UIND (J ,1 ) *360. 

102  CONTINUE 

1 REAO*,YLD,FF,GOL,NP,IEXP 
00  2 J=l,97 

00  2 K=4,6 

2 PR(J,K)=0.0 

TEST  FOR  USER  DIRECTIONS  I »=  CONTINUE,  0=  READ  NEW  YLO, 

-=  STOP 

IF(YLO)  99,131,3 

3 PRINT81,YL0,FF,G0L,NP,IEX» 

C7  COMPUTES  THE  MEAN  CLOUD  HEIGHT  0^  A SEA  LEVEL  PURST 
C2= 0.3048* ( 44. *6. 1* A LOG (YLO) -0.205* (A  LOG (YLO) *2. 42)* 

1 APS  (ALOG(YLO) *2.42) ) 

IF(C7.GT.20.)  GO  TO  5 

CRsl4.661*YL0**3.431 

SIGMA=1.*CR 

4 ACT=YLO*FF*2.035E*22 
GO  TO  8 

5 PRINT82,CZ 
IFINF.EQ.O)  GO  TO  1 
00  7 J=1,NP 

7 REAO*,PX,PY 
GO  TO  1 

8 00  10  J=l,97 

10  PR(J,3)=SIGHA* (l.»(J-  1)*0.020) 

PRINT83,YL0,FF,YLD*FF,G0L,CZ,CR,SIGMA,ACT 
COMPUTE  FALL  TIMES  AND  DISPLACEMENTS  FOR  EACH  GROUP 
I«CZ 
A«CZ-I 
A2«A/2.*I 

08*0.0  . 

15  00  2C  J*l,97 

CALL  FALL(A2,PR(U,1) ,A) 

PR(J,4)»PR(J,4)*nT 

PRU,5l*PR(J,5)*nx 
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20  PRC J,6)=PR{J,6)+0Y  ■- 

IFCOe.GT.O.a)  GO  TO  32 
«2«I«0.S 
00  3C  K=1,I 
«2*A2-1.0 

IF(A2.LT.<G0L*1.))  GO  TO  31  , 

00  3C  J=l,97 

CALL  FALL(A2,PR(J,1) , 1.) 

PR(J,4)  = PR(J,4  )*''T 
PR(J,5)=PP(J,5l*rx 

30  PR{J,6)=PR(J,6>»0T 
GO  TO  32 

31  OB=A2*0.5 
AsOB'GOL 

A2»(rB-G0L) /2.*G0L 
GO  TO  15 

32  PRIN’ee, < (ORCNjM) ,M=1,6) ,N=1,97) 

IF(NP.GT.O)  GO  TO  95 

C TEST  FOR  1500  KM  LIMIT 

33  00  48  J=l,97 

IF  (SORT  ( (PR  (J,  5)  -PR(1,5))**2KPR(J,5)  -PRC1,6))**21  .GT.ISOO.) 
1 GO  TO  49 
, 48  CONTINUE 

J»J-1 

'C  DEFINE  DISTANCE  LIMITS  OF  THE  PATTERN 
' 49  AHa<PR(J,6)-PR(l,6))/(PR(J,5)-PRtli5)) 

AN*5.*PR(J, 3) 

0X*(ABS(PR(J,5)-FR(1,5) )+5.*(PP{J,3)  ♦PPC 1 , 3) ) ) / 1 00 .+2. 0 
0V*(ABS(PR<J,6)-PR<1,6) ) *5  .♦ (PR( J, 3) ♦PR ( 1 , 3) ) > /1 00 . ♦2. 0 
OXL«(PR( Jt5> -5.*PRt J,3>> /OX 
OYL*(PR<J,6)-5.»PR(J,3) ) /OY 

IF(PP(J,5) .GT.PR(1,5) ) OXL=(PR(l»5) -PP (1 , 3) *5) /OX 
ir(PR(J,6) .GT.PR (1,6) ) OYL*(PR(l,5l-5.*PR(l,3>)/OY 
< 60  00  61  J-1,11 

j 61  OXX(J)=DX*(DXL»1+(J-1)»10.) 

PRINT84,CX,0Y, DXX 
ICKA=0 

BEGIN  OOSE  RATE  CALCULATTONS  FOR  EACH  POINT  (PX,PY)  OF  A 

rectangular  grid 

00  79  K=l,101 
PY»(0YL+131.-K)*0Y 
ICKB^O 

00  78  J=l,130 
PX«(nxL^J)*OX 
C THESE  6 CAROS  REOUCE  THE  CPA  TIME  REOUIREO  9Y 
C ELEMINATING  unnecessary  OOSE  RATE  CALCULATIONS 
YN0*AM*(PX-PR(1,5) ) ♦PR(1,6) 
AC*Sr'RT((PX-oR(l,5))*“2^(YN0-PR(l,6)  )**2) 

IFCAC.LT.O.l)  GO  TO  65 
AF»PY-YNO 

AE*APS ( (PX*PR (1, 9) ) ♦ AF/AC) 

IFCAE.GT.AN)  go  to  70 

65  CALL  FIELD (PX,PY) 

C OETER“INE  THE  OOSE  RATE  RANGE  AT  TM*  POINT  (PX.PY) 

IF(0R-1.)  70,71,66 

66  IF(0P-13.)  71,72,67 

67  XFCOP'IOO.)  72,73,68 
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68  IF<OR-1000.)  73,74,69 

69  IF(DP-3000.)  74,75,75 

70  rXY<J)=BA 
ICKB=ICKBH 
60  TO  78 

71  IXY<J)=B1 
ICK*=1 
60  TO  78 

72  IXY  (J)  =02 
60  TO  78 

73  IXY  (J)  =93 
60  TO  78 

74  IXYtJ)=B4 

60  TO  78  ( 

75  IXY<J)=95 

78  CONTINUE 

IF(CICKA.EQ.O) .AN0.<ICKB.GT.99n  30  TO  79 
IFCdCKA.GT.O)  .AN0.(ICK9.GT.99»)  30  TO  90 
PRINT85,PY,IXY 

79  CONTINUE 

THESE  STATEMENTS  .allow  coR  An  EXPA'^OEO  SCM.E  WHICH 
SHOWS  MORE  DETAIL  OF  THE  MAIN  PORTION  OF  T^JE  PATTERN 
' 90  IEXP»IEXP-1 

j IFdFXP)  1,1,93 

93  DX«Dy/2. 
j 0Y*0V/2. 

60  TO  to 

95  CALL  SINGLE (NP) 

IFCIEXP.GT.O)  GO  TO  33 
60  TO  1 

80  FORMAT (“1”5X“9ASE  DATA  INPUT~//76<">ARTI3LE  0ATA"25X“HIN0  DATA" 

1 //20(6X,F6.1,3X,Fi3.p,3X,F6.2,lOX,F5.1,3X,F5.1,3X,F7.5/) , 

2 77(6X,F&.1,3X,F10 .8,3X,F6.2/) ) 

81  FORMAT  (“l-SX^'INOUT  OATAi  “7F10. 3 ,216//) 

82  FORMAT  l//“  “5 <“X") BX^C^s  “F5.1”  <H,  ANO  EXCEEDS  20  KH.  PROCEED  " 

1 “TO  NEXT  PROBLEM"//) 

83  FORMAT <“  “5X~HEAP0N»S  TOTAL  YIELO  - “F7.3"  MT“//6X“FISSI0N  “ 

1 “FRACTION  - “F5.3/6X"MFAP0')*S  FISSION  YIELD  - “F7.3“  MT”/ 

2 6X“GRND  LYL  OF  THE  PATTERN  - "F7.2"  KM“/6X"MFAN  CLOUD  " 

3 “HEIGHT  - "F5.1“  KM“/6X“MEAN  CL  QUO  RADIUS  - "F5.1“  •CM"/6X 

4 “SIGMA  a "F5.2"  i<’M"/6X“T0Tai  GAMMA  ACTIVITY  AT  1 HOUR  = " 

5 1P£12.5“  GAMMAS  °ER  SECONO"//10 X ,40 ( “♦  ”)/10X“M0TE«  ALL  " 

6 "DISTANCES  ARE  IN  KILOMETERS,  AND  ALL  INTENSITIES  ARE  IN  " 

7 “R/HR  AT  1 HOUR“/10X,4a(”*  ")//10X"'1AP  LEGEND"//13X 

8 “1  - 1 TO  13  R/HR"/1CX“2  - 13  TO  103  R/HR“/10X"3  - 100  TO  “ 

9 “1000  R/HR"/10X"4  - 1000  TO  3300  R/HR"/t0X”5  - 3C00  OR  " 

A “MORE  R/HR“//) 

84  FORHAT(“l“9X"THe  HORIZONTAL  SCALE  IS  “F6.2“  TO  1 KM,  AND  THE  “ 

1 “VERTICAL  SCALE  IS  "F6.2“  10  1 KM." //10X“Tme  HORIZONTAL  “ 

2 “AXIS  REPRESENTS  THE  EAST-WEST  3IS?LACEMENT~//rx , 

3 10(F6.1,‘*X) ,F6. 1/1  ox, 10 (“•12345678 9”) “*“//) 

85  FORMAT (“  “F6.1,3X,100A1) 

86  FORMAT (“1"5X"PARTI:lE  GROUP  OATA“//5X“MIO-RANGE“5X“ACTIVITY“5X 

1 “STANDARD"9X"TIMF  T0“5X“EAST-WEST"5X~N0RTH-S0UTH”/6X 

2 “RADIUS  “5X“FRACTIcIn“5X“0EVIATI0N"4X“FALL"3X“DISPLACEMENT“ 

3 2X,“0ISPLACEMENT"/6X"(MICRaN9)"19X“(KH)"9X“(HRS.)“6X 

4 “IKH)“10X“(KM)“/e7(/6X,F7.2,6X,FlC.8,4X,F6.2,7X,F7.1,5X, 
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5 Pa.l,7X,F8.i)//) 

STOP 
END 


SUBROUTINE  FALL<H,P,OH) 

COHHON/BLOKl/HINn(?0,  T)  ,tJT,0X,0V 
DATA  PIR/l.TVSI’OESl 1PT9/ 

THIS  SUBROUTINE  lO^PUTES  FALL  TINES  AN3  OISPLACEHENTS 
FOR  FACH  GROUP 
U«l,A216E-05 

IFCH.LT.ll.)  U=1.789AE-05-3.3361E-07*H 

nn*h^i 

C0RS2.717866E- 13* WIND tNN , 31 *R**3' J» * 2 

Re*CPR/2A.-2.33&3E-04*COP»*2*2.0l5<»E-06»;OR**3-6.9105E-39*COR**4 
IF (OCR. GE. 138.)  R£=1C (-1. ERBTStO . 999* ALOGIO (COR) -0 . 046f»77* 

1 (ALOGIO (COR)) »»2»O.OC11235* (A  LOGIC (COR) )»*3) 

VT*1.8E*06*RE*U/(HINn(NN,3)*R)*(l.*0.1165/(R«WlNO(nN,3)) ) 
OT»OH/VT 

0X«-0T*MIND(NN,2)«C0S (WIN0(NN,1)  *PIR/100.) 

Oy*-DT*HINO(NN,2)*SIN (NINO (NN,1)*=IR/ 100.) 

RETURN 

ENO 


SUBROUTINE  FIELD (PX,PT) 

C0HMCN/9L0<2/PRtB7,6) , OR, ACT 
THIS  SUBROUTINE  CONPUTEB  AND  SUNS  THE  OOSE  RATE 
CONTPieuTEO  BY  EACH  PARTICLE  GROU» 

SUNxC.C 
OO  AOO  L=l,97 

EXe<(PX-PR(L,5))  •♦2*(P''-PR(L,6))**2)  /2./PR (L , 3) **2 
XF(EX.GT.;0.)  GO  TO  LCO 

SUM*SUM4PR(L,2)*C .3989*EXP(-EX)/»R(L, 3)**2 
400  CONTINUE 

OR=1.147125E-16*SUM*ACT 

RETURN 

ENO  . _ • . 


SUBROUTINE  SINGLF(N) 

COHHCN  /BL0K2/FR(97,6),0R,ACT 
|C  THIS  SUBROUTINE  CONFUTES  THE  OOSE  RATE  AT  A SINGLE  POINT 
I PRINT88 

i OO  600  NJsl.N 

REA0»,PX,PY 
CALL  FIELO(PX,PY) 

600  PRINT87,PX,PY,  OR 

$7  FORMAT (“0“5X"AT  THE  POINT  X*  “FS.l"  KN,  AND  Ys  ”F6.1~  KN,  THE  “ 

1 “RADIATION  INTENSITY  CUE  TO  '’A. LOUT  IS  “^7.1"  R/HR  AT  1 HR“) 

86  F0RMaTI“l“5X“INTENSITIES  FOR  SPECIFIC  POINTS  ARE  '.S  FOLLCHSI"/) 
RETURN 
ENO 
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